Bell inequalities with continuous angular variables 
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We consider bipartite quantum systems characterized by a continuous angular variable 9 £ [—7!", 7r[, 
representing, for instance, the position of a particle on a circle. We show how to reveal non-locality 
on this type of system using inequalities similar to CHSH ones, originally derived for bipartite spin 
1/2 like systems. Such inequalities involve correlated measurement of continuous angular functions 
and are equivalent to the continuous superposition of CHSH inequalities acting on bidimensional 
subspaces of the infinite dimensional Hilbert space. As an example, we discuss in detail one applica- 
tion of our results, and we derive inequalities based on orientation correlation measurements. The 
introduced Bell-type inequalities open the perspective of new and simpler experiments to test non 
locality on a variety of quantum systems described by continuous variables. 
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Introduction: The pioneering discussions of Einstein, 
Podolsky and Rosen (EPR) [1 , which rose the possibility 
of eventual conflict between the classical an the quantum 
definitions of realism and locality, dealt with measure- 
ments of continuous variables (CV) of a quantum system, 
as position q and the kinetic momentum p. Almost 30 
years elapsed before J. S. Bell promoted a regain of inter- 
est on the subject, deriving a quantitative criteria estab- 
lishing the border between quantum and classical physics 
concerning locality and realism [5]. The discretized ver- 
sion of the EPR paradox and of the Bell inequalities us- 
ing spin-1/2 like systems came with the work of Bohm [3 
and Clauser, Home, Shimony and Holt (CHSH) [?j and 
was motivated by the necessity to find simpler setups 
for experimentally testing the, up to then, Gedankenex- 
periments exposing the quantum-classical contradiction 
[5]. The almost 30 years following such pioneering ex- 
periments led to theoretical and experimental advances 
pointing out that using CV may provide advantages with 
respect to discrete systems. Higher intensity signals al- 
lowing to circumvent the detection loophole are among 
the announced advantages [6l [7] . Most of these works are 
based on canonical unbounded variables, as (q,p). They 
can correspond, for instance, to the sum and the differ- 
ence of two quadratures of the electromagnetic field. The 
first attempts to build CV non-locality tests consisted on 
finding ways to discretize the continuum to reuse con- 
cepts developed for the discrete case. It was shown in [8| 
that dichotomizing the phase space according to a state's 
parity and its displacement in phase space can lead to 
Bell type inequalities that can be violated by gaussian 
continuous variable entangled states. Other phase space 
dichotomizations are possible, as the one proposed in [9]. 
However, dichotomization does not correspond to a gen- 
uine CV measurement, since it provides an observable 
with a discrete, rather than genuinely continuous spec- 



trum. 

Reid and co-workers [7] proposed variance based Bell- 
type inequalities not demanding dichotomization, that 
are independent of the spectrum of the measured ob- 
servables. Another possible approach are entropy based 
inequalities that may be advantageous to detect entan- 
glement in non-gaussian states, since they go beyond the 
usual criteria that use second-order moments [10l [IT] . 

Up to now, most of the results on CV entanglement and 
non locality have been devoted to the case of canonical 
variables with an unbounded spectrum, as position and 
momentum. In the present work, we consider a different 
type of CV. We deal here with two quantum systems, A 
and B, characterized by angular variables 9i £ [— it, tt[, 
i = A,B on a circle, instead of a position q in the full 
real line. 9 can represent, for instance, the position of 
a particle on a circle, a rotator on a plane (for exam- 
ple, molecules confined to a plane [TJ]), or the polar an- 
gle locating the photon field in the plane transverse to 
its propagation, as illustrated in Fig. [I] The main re- 
sult of the present paper is to build inequalities involv- 
ing correlated measurement of continuous, real and pe- 
riodic angular functions fi{0i) (i — A,B), Surprisingly, 
the derived inequalities share similarities to CHSH-likc 
ones. As a matter of fact, their general form is given by 
a properly balanced continuous superposition of CHSH- 
type inequalities. Exploiting the possibility of revealing 
non locality using angular measurements is an approach 
that can be particularly useful and suitable for a num- 
ber of physical systems in atomic and molecular physics, 
photonics and mesoscopic systems in condensed matter 
physics. 

Theoretical model: As a starting point, we recall the 
CHSH inequalities for a bipartite two-level system. It 
can be written in a simple form as \{B(a, a 1 , b, b'))\ < 2, 
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with: 

B(a,a',b,b') = cr a <g)a b + cr a <g)a b > +a a > ®a b -a a > ®a b > (1) 

where o~ a is the Pauli matrix in one direction a of the 
tridimensional space, a and a' refer to different di- 
rections of the spin projection of subsystem A while b 
and b' are associated to subsystem B. The inequality 
\(B(a, a' ,b,b'))\ < 2 is fulfilled in the framework of lo- 
cal hidden variable (LHV) theories. The fact that there 
exist quantum states violating these inequalities disproof 
quantum mechanics as a local or EPR realistic theory. 
For quantum mechanical systems, ([!]) is bounded by the 
Tsirelson bound of 2y/2 [13]. 

We now consider a different situation of mea- 
surement providing continuous and 27r-periodic out- 
comes. In other words, the underlying Hilbert space 
is the square integrable 27r-periodic functions space 
(L 2 ({— 7r, 7r[, d9)). Such functions can be spanned by 
the basis {\mi};mi £ Z}, which are angular momentum 
eigenstates: J z \fni) = m\mi). An alternative continu- 
ous basis is {\9i);9i G [— 7r,7r[}, and it can be obtained 
from the \rrii) basis by Fourier transform: (rrii\9i) — 
jigimiOi^ j n ^ n j g re p resen tation, a local observable for 

each party can be written as Fi = J_ 7T d9ifi(9i)\9i)(9i\, 
with fiifii) real, bounded and periodic. 

Our guidelines to obtain continuous variables Bell in- 
equalities is to build a CHSH A\ Bell-operator similar to 
the one given by ([I]), but based on the correlated mea- 
surements of Fi for each particle. It is clear that, un- 
der the assumption that the spectra of Fi are bounded, 
this property is preserved by every unitary transforma- 
tion Ui(fa) such that F(fa) = Ui^^FiU} (fa). It is also 
straightforward to show that, for a LHV theory, we have: 

\{B{fa,<l>'ai<l>b,<l>b))\ = 

\{Fa{4>o) ® F B (fa) + F A {fa a ) <g> F B (fa) 

+F A (<f> a ) ®F B (fa b ) - F A (fa a ) ® F B {fa b ))\ < 2, (2) 

if the maximum value of fi(&i) is normalized to 1. What 
are the conditions the functions fi(9i) and the trans- 
formed operators Fi(fa) should satisfy to violate ([2| and 
allow for non-locality and entanglement tests? In order to 
answer this question, we focus on fi(Oi) = cos(#i) which 
is a function that can be measured on a variety of phys- 
ical systems. The results derived in the following can be 
straightforwardly generalized to all 2tt periodic functions 
such that f(9) = -f{9 - ir) V<9 e [0, tt]. 

Example of application: The observables C — Fi cor- 
responding to fi(9i) = cos(0j), can be expressed on the 
|m) basis as: 

C =J2 l(\m + l)(m\ + \m)(m+l\). (3) 

For a particle rotating on a circle, this operator corre- 
sponds to the projection on the x axis of the particle 



position; it can also be related to the spatial orientation 
of a two dimensional rotor, or to the phase of a supercon- 
ducting circuit. Operator C has a spectrum of doubly de- 
generated eigenstates \9) and |— 9), both with eigenvalue 
cos 9. For this operator, we can define an equivalent to 
the rotation of a spin 1/2 system: it is the unitary op- 
erator e lJ ^, (that is, the free evolution operator during 
a time t = 2%fa. for a free particle with unit mass and 
angular momentum J z , constrained to move on an unit 
circle). From C, we can define C(fa) = Fi(4>i) as: 

C{fa) = e lJ ^Ce- lJ ^ = \Yj e i(2m+1) + 1| + 

e- l(2m+1)< V + l)(m|, (4) 

and for a bipartite system (particles A and B) the 
Bell operator (p]) reads: B(<f> a ,fa a ,fa,(i) b ) = C A (<p a ) <£> 
C B (fa) + C A (faj® C B (fa) + C A ® C B (fa b )- C A (fa a ) ® 
C B (fa b ), where C A r B \(fa) are operators defined as in 
Eq. Q acting on the Hilbert space of particle A(B). The 
spectrum of C(fa) does not depend on fa Diagonalizing 
C{<t>) shows that its spectrum is bounded, with |(C)| < 1. 
Thus, for a LHV theory, \{B(fa,fa a ,fa,(j)' b ))\ < 2 holds. 
However, for some set of phases fas and some quantum 
states, this inequality is violated. In order to show that, 
we calculate the spectra of B(4* a , fa a , fa, fa b ) operators. 
Actually, such spectra depend only on the relative phases 
fa - fa a and fa - fa b HU US]. We can define B(£ a , &) = 
B(fa a -fa,fa b -fa) = B(0, fa a -fa,0,fa b -fa) = e lJ ^ ® 
e lj2 B^B(fa,fa a ,fa,fa b )e- u A^®e- u B^K B(£ ,£ b ) and 
B((j> a , fa a , fa, fa b ) are related by an unitary transforma- 
tion. Therefore, the variation of only the two phases 
(£ai£b) is enough to explore the spectrum of all the 

B(fa,,fa a ,fa,fab)- 

As a first approach to the study of the dependence 
of operators B(£ a ,£ b ) with the phases £j , we discretize 
the possible outcomes of the correlation measurements. 
For this purpose we consider that the measured quantum 
states lie on a 2M+1 finite dimension space generated by 
the basis set {|m); m G Z, \m\ < M}. We thus define the 
projection C( M ^(£) of the C(£) operator on this reduced 
space as : 

CW(£) = i ^2 e l(2m+1) t\m){m+l\+e- l[2m+1 ^\m+l){m\ 

\m\<M 

(5) 

and C (M) = C (M) (0). We also define discretized Bell 
operators Z?( M ) and B^ M \ analogously to the continu- 
ous case, but where the C(£i) operators are replaced by 
C^ M \£,i) operators, which act on the considered finite di- 
mensional space only. Discretization allows to explicitly 
compute the 2M+1 eigenstates and corresponding eigen- 
values of C< M \ providing thus a more intuitive physi- 
cal image of the considered operators. Also, it enables 
the numerical study of the spectrum of £?( M ) (£ a , £ b ) and 
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FIG. 1. Examples of physical systems for which our results 
can be applied: (a) The transverse profile of a propagating 
light beam (b) Particles moving on a circle (as charged parti- 
cles in an electromagnetic trap) (c) rigid rotators confined to 
two dimensional motion. 




FIG. 2. Plot of the maximum eig envalueof£ (M) (£a,6),& 

max , 

in the region of violation > 2) as a function of £ a and 

£b for M = 2 (inner blue plot), and M = 5 (outer red plot). 



its phase dependency. Denning &^ x (£ a ,£a) as the high- 
est eigenvalue of B^ M \^ ai ^b) we see in Fig. [2] that, for 
M = 2 and M = 5, &max(£ a ,£a) reaches its maximum 
at £ a = £(, = tt/2. We verified numerically that this 
fact is independent of M for 1 < M < 20. In addi- 
tion, |6(7r/2, 7r/2)| increases as M increases, being already 
greater than 2 (i.e., enabling the Bell-type inequality vi- 
olation) for M = 2. 

The study of the discretized operators is now used as a 
starting point to the derivation of the main results of this 
Letter, which are non-locality and entanglement tests in 
the continuous limit. We start by analyzing the point 
£ a = £,b = t/2, since it is the one that provides the 
highest contrast between 6 m ax(£a,£a) and the violation 
threshold value (2). 

We then move to the continuous limit (M — > 00) and 
search for the highest possible value of violation and the 
corresponding non local states for this choice of phases. 
Using the expression of C(£) Q and the anti-symmetric 



property of cos0, we compute (0'|C(f )|0) explicitly: 

/>7T 

C(-)= / d0cos0[i|0-7r)(0|-i|0)(0-7r|]. (6) 
2 Jo 

It is now convenient to define the following operators : 

(7) 



a*(9) =|0)(0| + |0)(0| 
°»(0)=i(\e){9\-\0){0\) 
a*(9) = \9){9\-\9){9\, 



for £ [0, 7r[ and with = — ir. These operators consti- 
tute an orthogonal continuous set of Pauli-like operators, 
and they fulfill relations analog to the usual Pauli matri- 
ces: ai(9)<7 k (0') = ia l (0)6{9 - 0') where {j,k,l) denotes 
circular permutations of (x, y, z). They also fulfill orthog- 



onality relations a j (9)a j (6') = S{9 - 0')(|0) 
Eq.(|6]) can thus be recast as: 



C 



d0COS0CT Z (0) 



C(j> = 

Thus, the Bell operator B n 
the following direct sum: 



d0 cos 9a y {9). 



?> (op- 
ts) 

(9) 



§ ) can be written as 
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d0d0' cos cos 9'X(9,9') 



(10) 



where 



X(0, 9') = [a A (9) ® aW) ~ <t v a V) ® °W) 
-a A (9)®<(0')~°\(0)®<(0')], (11) 

and 0^(0) and <J 3 B (9) (j = x,y,z) are operators defined 
by Eqs. ^ on the Hilbert space of parties A and B re- 
spectively. The X(9, 9') operators are orthogonal: 



X(9 a , 9 b )X(9' a , 9' b ) = 4<5(0 a - 9' a )S(9 b - 9' b )x 

(\9 a ) (9a\ + \9~a) (~9~a\) ® (|^> (9 b \ + \¥ b ) (9~ b \) 



(12) 



and completely analog to the usual 2-qubits (4- 
dimensional) CHSH operators. We have thus the sur- 
prising result that the Bell operator B m is the weighted 
continuous direct sum of 2-qubit-like CHSH Bell opera- 
tors X(9, 0'), with the weights being given by cos cos 9' . 



Thanks to the orthogonality property (12 1, finding the 



spectrum and eigenstates of B m is a simple task. Indeed, 
for each 9 and 0' it is enough to diagonalize the 4x4 
matrix representing X(9,9'). We find that B m can thus 
be written in diagonal form as: 



B„ 



2V2E 



n=±l 



d0d0' X 



cos0cos0'| X "(0,0')}(x'W)l (13) 
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where 



\x tl (W)) = j^l\6)®\(>') + \0)®\0') 

Ti(V2Ti)(\e)®\W) + \e)®\e'))} 



(14) 



are the eigenvectors of X(8, 6') and of B m with non-zero 

1 /2 

eigenvalues and where N± = [2(2 =p v2)l is a nor- 
malization factor such that : ( X n ' (6' a , 0' b )\x n (0 a , 6 b )) = 

Therefore, the spectrum of B m is continuous and equal 
to [-2\/2, 2 \/2] • I (-Bm) I is thus bounded by 2^/2 as in the 
2-qubit CHSH inequality. However, the values ±2y/2 are 
only reached by eigenstates of -B m which are not physi- 
cal states since they are composed of perfectly oriented 
states, namely |0) and \tt). In the angular representa- 
tion, these states are expressed by delta functions, since 
(8\#) =5(6-6'). 

In order to explicitly construct physically sound states 
violating our Bell inequality, we should consider wave 
packets consisting of continuous superposition of eigen- 
states \x +1 (6,9')), with 9 and 9' localized around 6 — 
9' = 0, point of the maximum eigenvalue of B m . An 
example of such wave packet is given by: 



I*) 



d9 / d6'g a (8)g b (6')\ X +1 (6,6')). (15) 
Jo 



where g a (6) and gb(9) are normalized L 2 (]0, ir[, dO) func- 
tions with support containing 9 — 0. The expectation 
value of B m for this state is 



(Br, 



2V2 



d6 



d6' cos 8 cos 9' \g a (9)\ 2 \g b {9')\ 



/\|2 



(16) 

The wave packet ( 15 1 can be produced making a linear 

(17) 



combination of the one-particle wave packets 



|S) = / d6g(6)\8) and \g) = / d6g(6)\8), 
Jo Jo 



taking the same coefficients of Eq. ( 14 ) 
1 



I*) = ^f^a) ® \9b) + \ga) ® |^» 

- i(V2 - 1) (\g a ) ® \gj) + \g$ ® \g b ))} 



(18) 



It is useful now to obtain a simple relation between the 
wave packet localization and the violation of the Bell 
inequality. For such, we take for g(6) the ideal case of an 
angular slit with aperture 69, given by : 



9a(0)=9b(0) 



-4= for 6 < 69 

Jse 19 

otherwise 



Using Eq. (16 1, the value of (B m }q, can be written as 
function of the aperture 69 of the slit as follows: 



(B m ) 9 = 2V2 



sin < 



56 



(20) 



This equation above shows that we can obtain values of 
(B m )\B > 2, violating the Bell inequality, with an aper- 
ture of 69 < 18° which is not a too restrictive condition. 
There is thus a relatively broad collection of simple two- 
particles non local pure states involving coherent super- 
position of one-particle wave packets localized around 
9 = and 9 = ir that violate the derived Bell inequality. 

We now study the example of even more real- 
istic states, which are non-pure ones, establishing 
some conditions for them to violate the derived Bell- 
type inequalities. We consider the analog of the 
Werner states [To] : p(r}) = r\pA ® Pb + (1 — 



where p A ( B ) 
+■ 1 0(1(6)) (da(b) 



= Tr 



B(A) 



l*>(*|] = 



|), with I*) given by ([18 1. 
6b)) does not depend on 



77)|*) 

\ (|fa(b)) (g a (b)\ 

As (x n (9a,8b)\PA®p B \x n 
n = ±1 for all 6 a ,6 b G [0, n], then Tr [B m p A ® p B ] = 0. 
It turns out that the expectation of B m for the state 
p(rj) is simply given by (1 — rj)(B m }y. In the ideal case 



where g a and gb are given by (191, the maximal allowed 



value of the mixing coefficient 77 for the Bell equation 
to be violated is a simple function of the slit aperture : 



Tr [B m p{rj)\ > 2 => r) < 1 



1 

V2 



se 

sin S8 



) , providing condi- 



tions for this type of mixed states to violate non-locality. 
Conclusion: we derived general CHSH-like Bell-type in- 
equalities for continuous variables using measurements of 
bounded observables with anti-symmetric and periodic 
spectra. Such observables are relevant for a number of 
experimental systems, ranging from the photonic orbital 
angular momentum to the movement of material parti- 
cles on a circle and the phase of superconducting cur- 
rents. We discussed in detail an example of observable, 
but the results obtained are valid for all observables with 
the same spectral properties and symmetries. Our results 
exploit angular variables from an original approach, and 
open the path to novel experiments with a wide range 
of quantum devices, throwing some light on the difficult 
task of non-locality violation with genuinely continuous 
variables. 
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